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■ Abstract 



In a recent paper (Phys. Rev. A78, 020101(R) (2008)), Kim et al. have reported a large anomaly 
in the scaling law of the electrostatic interaction between a sphere and a plate, which was observed 
Qh| during the calibration of their Casimir force set-up. Here we experimentally demonstrate that this 

C ' behavior is not universal. Electrostatic calibrations obtained with our set-up follow the scaling 



^3 . 



law expected from elementary electrostatic arguments, even when the electrostatic voltage that 
one must apply to minimize the force (typically ascribed to contact potentials) depends on the 



' separation between the surfaces. 
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Casimir force experiments are routinely used to set new limits on Yukawa corrections to 
the Newtonian gravitational attraction between surfaces at submicron separation (see |l| 
and references therein). To strengthen those constraints, new experiments must provide 
data with errors smaller than any other previous measurement. With claims reaching 0.19% 
relative experimental errors at a 95% confidence level it is important to ask whether 
there exists any technical challenge or physical mechanism that might impede any further 
improvement in this important field. 

In a recent paper Kim et al. have reported a systematic effect observed during the 
calibration of their Casimir force set-up that might represent a severe limitation to the de- 
velopment of future experiments on the Casimir effect, and that, in some extent, might have 
been overlooked even in previously reported high accuracy measurements. They observed 
that the electrostatic force gradient between a ~ 30 mm radius spherical mirror and a metal- 
lic plate scales with surface separation d like ~ which represents a 15% deviation on 
the exponent with respect to the l/d? behavior expected from elementary electrostatic cal- 
culations. If this anomaly were confirmed to be a general phenomenon for metallic surfaces 
at very close separations, all the arguments used to calibrate high accuracy Casimir force 
set-ups, which relies on elementary analysis of the electrostatic attraction between the two 
surfaces, would be invalidated with severe consequences on the results of those experi- 
ments. Driven by these considerations, we have performed a high precision experiment to 
investigate the electrostatic force between a sphere and a plate in the separation range from 
~ 100 nm to ^ 2 /im. 

The main goal of this paper is to test the validity of: 

^_ SonRiV + Vof^ 



d 

where F is the electrostatic force expected between a spherical surface of radius R and a 
plate kept at a separation d with d << R, and where Eq is the permittivity of vacuum, V is 
the applied voltage, and —Vq represents the voltage that one has to apply to obtain minimal 
electrostatic force (typically ascribed to contact potentials) 

In the apparatus used for this experiment (see Fig. [T]), the sphere is directly glued under 
the hanging end of a micromachined cantilever (the force sensor), which is then mounted 
inside the measuring head of a commercial atomic force microscope (AFM). The plate is 
anchored to a custom-designed mechanical stage, which is fixed underneath the AFM head. 



The mechanical stage allows one to bring the plate in close proximity with the sphere, and 
to perform measurements of the force as a function of separation. 

The force sensor used in this experiment is a rectangular Si cantilever (525 x 35 x 4 fim^) 
with nominal spring constant k ~ 0.9 N/m. A 100 fim radius polystyrene divinylbenzene 
sphere (Duke Scientific) is attached to the free end of the cantilever with UV curable glue. 
The cantilever and the sphere are then coated with a 10 nm Ti adhesion layer and a 200 nm 
Au film by magnetron sputtering. A similar coating is also deposited onto the plate - 
a 5 X 2 mm^ polished sapphire slide. The AFM head is a low-noise Veeco Multimode, 
which exploits standard optical lever techniques to measure deflections of the cantilever 
with 0.1 nm precision over a 50 kHz bandwidth. The mechanical stage consists of a stick- 
slip motor (Attocube) and a piezoelectric translator (PI). The stick-slip motor is moved only 
at the beginning of the experiment to bring the plate within a few microns from the sphere, 
while the separation between the two surfaces during the actual experiment is varied with 
the piezoelectric translator. This translator, which is controlled by a capacitive feedback 
loop, has been calibrated by the manufacturer against traceable standards. Its resolution is 
reported to be equal to 50 pm. The mechanical stage and the AFM head are anchored to a 
1 dm^ aluminum block that is maintained at a temperature ^ 10 K above room temperature 
by means of a feedback controlled heating system. The block is mounted on a commercial 
active anti- vibration stage (Halcyonics), which is placed inside a 1 m^ acoustic isolation box. 
Finally, this box is placed on an optical table inside a temperature controlled room. 

In order to measure the dependence of the electrostatic force between the sphere and the 
plate as a function of their separation, we apply an oscillating voltage V = Vdc+Vac sin (ut). 
Using Hooke's law and Eq. [H the signal S of the photodetector of the AFM head can be 
written as 



^ = 7' 



k 



kdP 



(Vo + Vncf + ^ + 2 (Vo + Voc) Vac sin {ut) - ^Ac^os{^-^t) 



(2) 



where k is the spring constant of the cantilever and 7 is a parameter that translates cantilever 
deflection into signal (expressed in V/m). Note that we have added an exponent p to the 
denominator: our goal is to verify whether p = 1 or not. The signal S contains two static 
components and two periodic components at angular frequencies uj and 2uj. The static 
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components of the signal represent the static deflection of the cantilever, which, as we will 
show later, is always smaller than 0.2 nm, and will thus be neglected. The two periodic 
components can be used to fully characterize the electrostatic interaction. For this reason, 
the photodetector is connected to two lock-in amplifiers operating at frequencies uj (LI) and 
2uj (L2), as illustrated in Fig.[T] From Eq. [21 one can see that the output of LI is proportional 
to Vo + Vdc- One can thus create a negative feedback loop where LI generates Vdc in such 
a way that it keeps Vq + Vdc small, as typically done in Kelvin probe force microscopy j^. 
The loop gain G of the current experiment varies from 10'^ (at ~ 2/im separation) to 10*^ (at 
~ lOOnm separation). Because Vdc = ~(^Vb and \ Vdc\ < 50mV in all our measurements, 
the feedback loop certainly compensates Vq down to \Vdc + K)| < 50 fiV. Thus, one can 
assume Vdc = —Vq for all practical purposes. We stress that the purpose of this feedback 
loop is two-fold: the compensation voltage Vdc is measured accurately at all distances and 
the static defiection of the cantilever is minimized by effectively zeroing the first term of the 
expansion in Eq. O As far as L2 is concerned, note that its peak-to-peak value is given by: 

where a is proportional to the curvature of the parabola described by Eq. [H which can be 
obtained as S2uj/V^q. Therefore, by examining the measured values of a as a function of d 
one can verify whether p = 1 or not. To obtain a as a function of d, we start by placing 
the plate a few micron away from the sphere. We then move the plate towards the sphere 
in discrete steps. For each position, we measure 820; for a properly chosen value of Vac- At 
first, one might think to simply keep Vac constant during the whole run. However, it is more 
convenient to reduce Vac such that 82^; stays constant as the surfaces approach, because, in 
this way, the relative random error of a is equal for all separations: aa/d = (^32^/821^ (if we 
assume that the random error on Vac is negligibly small). Therefore, before we move the 
plate to the next measurement point, we use the values of a measured in the same run at 
larger separations to estimate the value that a should assume in the next position, and we 
set the new value of Vac accordingly 61] . This procedure does not by any means introduce 
systematic errors on a, which is still evaluated as a = 82uj/Vac, where 82^ is the actual value 
measured by L2. Systematic errors on Vac can also be ruled out, because we calibrated the 
digitally controlled function generator before starting the experiment. 

In the current experiment, L2 has a 24dB/octave low pass filter with a 1 s RC time, 



and data are acquired with 5 s integration time at every position. We use an 82^ set-point 
that corresponds to a peak-to-peak cantilever movement at 2u; of roughly 0.3 nm [3]. With 
this set-point, Vac varies between 450 mV at (i ~ 2 /im to 100 mV at (i ~ 100 nm. This 
corresponds to a static deflection of the cantilever at Vdc = — Vq of roughly 0.15 nm, which 
is thus negligible. To minimize drifts of the amplitude response S2U) due to potential changes 
of the resonance frequency of the cantilever, we work in the quasi-static regime, setting uj/2% 

resonance frequency of the force sensor (1.65 
kHz, as obtained with an independent measurement). The total measurement consists of 
184 runs over 1050 minutes. 

In Fig. [2] we show the value of a as a function of the extension of the piezoelectric stage 
dpz (see Fig. [T]) for one randomly chosen run. If we neglect the static deflection of the 
cantilever, the actual separation between the surfaces d is given hy d = — dpz, where do 
is the initial separation, which is a priori unknown. To validate whether p = 1, we fit the 
data with an equation of the form 

[do — dpzY^ ^ ^ 

where do and k = f^ee parameters. The fits are performed using standard 

minimization algorithms, for which it is necessary to first estimate the error on the data. 
This is done by measuring a at a fixed dpz for 120 minutes. The results are shown in 
Fig. [3l The data distribute along a smooth curve that is not constant because of drifts in d^ 
and/or k |9j]. In the inset, we plot a histogram of the relative difference between the data 
points and the smoothed curve. One can clearly see that a follows a normal distribution 
with a standard deviation of 0.56%; the relative error in a single measurement of a is thus 
0.56%. Note that this error represents an uncertainty of 600 fm in the determination of the 
root-mean-square motion in response to the varying potential difference Vac %■ Using this 
.elat™ e.o, we repeat the Bt thtee t.„tes: lettmg p a= a ttee parameter, and forcng p = 1 
(as elementary electrostatic calculations suggest) or p = 0.7 (as found in [Z\). If p is a free 
parameter, one obtains p = 1.005 ± 0.004 with reduced equal to 1.19. The fit with p = 1 
gives rise to a comparable value of reduced ix^ = 1.21). The fit with p = 0.7, on the 
contrary, produces a reduced of 411. Our data thus follow the behavior expected from 
elementary electrostatic arguments. It is important to stress that these values are obtained 
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analyzing all data except those for which d < 120 nm (see the open squares of Fig. [2]). If 
those data are included, the quickly increases. This is not surprising because, at small 
separations (smaller than ~ 120 nm), the Casimir force bends the cantilever so much that, 
within the precision of the current experiment, one cannot assume d = — dp^ anymore. 

To make our claim more robust, we fit each single data set with p = 1 and analyze the 
behavior of the reduced x^- For the sake of computational convenience, we rewrite Eq. [Has 

1 1 

- = - {do- dp^) . (5) 
a K 

To avoid systematic errors due to the bending of the cantilever at small separations, we 
apply again a mask to the closest A^ — 41 data points, where N is the total number of points 
in a single run (similar to Fig. [2]). The relative error on a is so small, that 1/a also follows a 
normal distribution (with 0.56% relative error). The average reduced over 182 runs (two 
runs are outliers with values of 4.8 and 2.9) is 1.03, and the values are distributed 
with a standard deviation of 0.23 ^]. The anomalous scaling law observed in is thus 
not a universal behavior. 

Now that the p = 0.7 behavior is ruled out, it is interesting to plot the residuals of the 
fits with p = 1. From this plot, reported in Fig. HI it is evident that our data systematically 
deviate from the fit, with maximum deviation of ^ 1%. This behavior is due to the fact that 
eq. [1] is based on the use of the so-called proximity force approximation (PFA) [ll| . Using 
;he complete analytical equation for the electromagnetic force between a sphere and a plate 
1^ to calculate the average residual expected from a fit with eq. [5] (with k and d^ as free 
parameters), one obtains the continuous line of Fig. HI Note that, although the difference 
between the complete analytical equation and the PFA goes to zero as separation decreases. 
Fig. m does not show this behavior. This is due to the fact that we are fitting a linear 
function (eq. [5]) to slightly curved data. From the same figure, one can also see that the 
random errors at larger separations are actually smaller than that used in the fits (0.56%), 
which was measured at ~ 150 nm (see dashed line of Fig. H]). This justifies the fact that 
the fit with eq. [1] and p = 1 gave ~ 1 even if a not completely correct theoretical model 
was used. 

We want to stress that such an excellent agreement can only be obtained if the experi- 
mental apparatus is exceptionally stable. In the inset of Fig. HJ we plot the initial separation 
do obtained from the fits of all 184 runs. The total drift in do is 45 nm over 1050 minutes. 
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This means that our set-up drifts only ~ 40 pm/min (compared to 1 nm/min in jis!]), or 
0.24 nm drift per measurement run (compared to ~ 60 nm per run in |2|). From the fits, we 
also observe that the other fit parameter k drifts 1.1% over 1050 minutes. This corresponds 
to 0.001%/min or 0.006% per run, and is likely due to a change of 7 caused by a slow drift 
of the laser spot over the photodetector of the AFM head. The effects of both drifts are 
negligible in one measurement run. 

The authors of [2I noticed that, in their set-up, the voltage needed to minimize the force 
depends on the separation between the surfaces. In Fig. [5] we plot Vdc as a function of d for 
the data set shown in Fig. [21 where was determined by fitting a as a function of dp^ with 
p = 1. Also in our measurements, the compensation voltage clearly depends on d, varying 
by ^ 6 mV over 2 fim. This behavior is reproduced in the other 183 measurement runs. In 
order to show that the dependence observed is not due to drifts of Vdc with time, in the 
inset of Figj2] we plot the value of Vdc at (i ~ 275 nm as a function of run number. It is 
clear that, at this separation, Vdc drifts, at most, 40 fiV per run. A similar behavior is 
observed for all separations. 

Interestingly, the data seem to distribute along a curve that goes like a log d + b (reduced 

= 0.8). A similar behavior has been recently reported in an experiment between Ge 
surfaces [1^]. A rigorous explanation of the origin of the dependence of Vdc on d goes 
beyond the purpose of this paper. Still, our measurements, together with the results of |l^ 
and [2I], suggest that it is indeed of fundamental importance to measure Vdc for all values 
of d, as in 2l], and, previously, in 151]. Furthermore, since our calibration follows elementary 
electrostatic arguments, we conclude that, in general, it is not sufficient to check the scaling 
of a with d to rule out the presence of a distance-dependent Vdc- 
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FIG. 1: Schematic view of the experimental apparatus. 
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FIG. 2: Plot of Q (see Eq. [3]) as a function of the position of the piezoelectric translator for 
run number 107. The error bars on the data are within the dimension of the symbol. Black 
squares indicate data that are used for the analysis. Open squares are data that are excluded from 
the analysis. The continuous line shows the fitting curve obtained with elementary electrostatic 
arguments (a oc 1 /d). The dashed line represents the best fit obtained on the basis of the anomalous 
behavior observed in [2] (a oc l/d?''^). 
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FIG. 3: Plot of a (see Eq. [3]) as a function of time for a fixed value of d^z- The thick line represents 
a smoothed curve that accounts for the drifts in our set-up during this measurement. The inset 
shows a histogram of the relative difference between the data points and the smoothed curve, 
together with a Gaussian fit. 
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FIG. 4: Residuals of the fits with p = 1 plotted as a function of separation. The continuous line 
represents the expected deviations due to the use of the proximity force approximation instead of 
the whole analytical equation. The dashed line indicates the separation at which the data of Fig. 
[3] were taken. Inset: value of the initial separation between the sphere and the plate as a function 
of run number as obtained by fitting the data on the basis of elementary electrostatic arguments. 
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FIG. 5: Electrostatic compensation voltage as a function of sphere-to-plate separation for run 
107. The continuous line represents the best fit with Vdc = alogd + b (reduced = 0-8, 
a = —4.4 lb 0.2 mV, b = 4.3 it 0.6 mV, d in nm). The initial separation between the sphere and 
the plate is obtained from the continuous line of Fig. [5J Each error bar represents the standard 
deviation of the gaussian distribution of the 184 values of Vac at that separation. Inset: Vdc at 
d ~ 275 nm as a function of run number. 
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